Interlayer Exchange Coupling Beyond the Proximity Force Approximation 
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Ion bombardment has been shown to be capable of enhancing the interlayer exchange coupling in 
a trilayer system that exhibits giant magnetoresistance. We demonstrate that this phenomenon can 
be derived from the phase coherence among scattered paths within the two rough interfaces when 
their topographies are correlated. In the case of mild corrugations, our method reproduces the 
predictions by the proximity force approximation which does not consider the interference. When 
the characteristic Fourier conjugate of the tomography becomes large and comparable to the Fermi 
momentum, interesting new features arise and can only be captured by our more general approach. 
Among our findings, the scenario of an enhanced interlayer exchange coupling due to the interface 
roughness is explained, along with how it depends on the sample parameters. An additional channel 
for the resonant transmission is identified due to extra scattering paths from the roughness. 

PACS numbers: 43.30.Hw, 75.70.-i, 68.35.Ct, 79.60.Jv 



I. INTRODUCTION 

Interlayer exchange couphng (IXC) has been studied 
for more than twenty years [l|-|6[ with applications in 
phenomena such as the Giant Magnetoresistance [3, [^ 
(GMR) and the Tunneling Magnetoresistance (TMR). 
Due to the lack of reliable microscopic theories, the in- 
terface roughness (IR) was mostly treated by the static 
average [3, m, llfl which invariably led to a suppression 
on IXC. Improvement has been achieved by a systematic 
study using the perturbation method [11|. We shall follow 
up this line of approach with more detailed calculations 
and report new information on how to raise the sensi- 
tivity of GMR and why an enhancement in coupling is 
possible, as has already been observed in ion-bombarded 
samples |12l|. 

Besides IXC, the Casimir effect between metallic 
mirrors [13 faces the same complexity because a lot of 
experiments were performed by using a spherical test 
body in addition to the unavoidable corrugations on its 
surface. This Casimir force has been measured with 
high precision which provides a fertile ground to test 
the theoretical models of IR. The proximity-force ap- 
proximation (PFA), equivalent to statically averaging 
over the plane-plane geometry, is the first intuitive the- 
ory to be examined by both experiments \\A. [l5l | and 
theories fig. Since the PFA was shown to work only 
for mild corrugations |la| . Maia Neto et al. generalized 
it by the perturbation theorv flTHigll t o obtain better 
agreements with the experiments [20, [2l[. Recently, a se- 
ries of experiments on severe corrugated mirrors, which 
are outside of the applicability regime of the PFA, con- 
firm that a scattering approach such as Maia Neto's is 
needed to capture the essence of the nontrivial diffrac- 
tion effects. [11 [H] 

We observe that IXC shares the same physical con- 
cept and mathematical construction as the Casimir ef- 
fect, which connection is proved in Appendix A. A quick 
way to convince oneself of this similarity is by the expres- 
sion of the Casimir energy [2J, [23 E between two parallel 



fiat mirrors separated by a distance D and with area A 
and reflection coefficient r(fc|| , E) 
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where ^ is the imaginary frequency and p denotes the 
transverse electric and magnetic modes. This bears great 
resemblance to the IXC energy in a trilayer system with 
two metallic layers and a metallic or insulating spacer [j| 
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where fc|| integration is over the interface Brillouin zone 
while the D now represents the spacer width. One pur- 
pose of this work is to generalize this comparison to in- 
clude the effects of IR. However, the physics in trilayers 
is more versatile because the spacer exhibits a character- 
istic length scale, the Fermi wavelength 1/kp, compared 
to which other length scales can be tuned to give differ- 
ent behavior. These parameters include the roughness 
amplitude, corrugation length, and spacer width. Fur- 
thermore, the fact that the spacer can be either metallic 
or insulating also enriches the phenomenon caused by IR. 

We apply the perturbation method to calculate IXC 
in Section II, and demonstrate that the predictions are 
equivalent to those by PFA in the limit of smooth cor- 
rugations in Section III. In Section IV, the two interface 
tomographies in TMR and GMR are assumed to corre- 
lated, with special attention to the interference effect on 
IXC. Section V is devoted to study why, when and how 
much the IXC can be enhanced by IR. Discussions and 
conclusions are arranged in the final Section VI, where 
improvements over our previous work are explained. To 
preserve the conciseness of the main text, a rigorous proof 
of the connection between IXC and the Casimir effect is 
arranged in Appendix A. 



II. PERTURBATION METHOD 

In this section, we start by deriving the two- 
dimensional scattering states generated by a right- 
moving plane wave with momentum {kx,ky) that inter- 
acts with an irregular interface A{ij) at x = 0. The po- 
tential on the left side of A is set to be higher in energy 
by Vf). The wavefunctions on the left and right sides are 
denoted by $(x, y) and ^{x, y), respectively. The bound- 
ary conditions are 
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In this work we shall assume that A is much smaller 
than both l/k^ and the major Fourier corrugation wave- 
length Ac in order to proceed with the perturbative calcu- 
lations as in the Casimir effect |18j . The scattering states 
can be obtained by treating the IR as a perturbation to 
the smooth interfaces, 
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where ^oix,y) and 'I'o(a;,y) are the unperturbed scat- 
tering states, and the transmitted wave \I'o(a;,2/) carries 
momentum (fc^, ky) and superscript (n) denotes the n-th 
order perturbation. For an elastic scattering, the disper- 
sion relation in Eq.Q is 
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where to* represents the effective mass of the carrier. 
Insert Eq.(|l]) into Eq.© and use |fca;A(y)| and \k'^A{y)\ 
as the perturbation factors to expand Eg . (Bl) . Retaining 
up to the second order, one can show that|ll| 
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where the subscripts ky , qy denotes scatterings from ky 
to qy state, and tk^,k' is the transmission coefficient for 
a smooth interface. 

Interlayer exchange coupling in a trilayer system, af- 
fected by the quantum interference among the reflected 
waves, can be described by the reflection matrices Q 
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where W is the length of the interface, / is the unit 
matrix, RJ^ /R]^~ are the reflection matrices from the 
left/right smooth interfaces. 

When the topography AL{y)/Afi{y) at the left /right 
interface is considered, the reflection matrix can be writ- 
ten in powers of the perturbation: 
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The zero-order matrix R 

interface and is diagonal in the basis: 



corresponds to a smooth 
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where rr^.i. y lrj.„ „> are the reflection coefficients of 
momenta kxjqx from the left interface while k'j^^/q'j^^ de- 
notes the momentum in the left side layer. Same for the 
definition of R^- 

The first and second-order matrices is constructed by 
Eq.© and Eq.® as 
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Inserting Eq. (fT2)) into Eq.®, we can compute IXC up to 
the second order in A^ and Af(: 
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The first-order energy correction will be zero because it is 
proportional to the averages {Al) and {Aji) which are set 
to be zero by construction. The major correction, there- 
fore, comes from the second order perturbation and can 
be separated into correlation and uncorrelation terms: 
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Note that these results are general, which will reduce to 
Eq.(3) in Ref.[ll| when restricted to the double limits 
of kp\c S> 1 and kpD ^> 1, i.e., a wide spacer. The 
reason is that kpD 3> 1 allows us to ignore the energy 
contribution from higher order round-trip reflections in 
Eq.dlSl) and Eq.([TCl): 



5l\E, 



(2) 



W 



„^ f'^" dE f dky V-^, , 

qy 



'^ \^Rx ^ ^x)iL;k^,k'j,^iR\q^,qR^e. 

X {ky\AL\qy){qy\AR\ky) (18) 



5AE, 
W 



(2) 



27r 



2Im ^ 

{fc;,(A|)-^(,^ 



rf£' /■ dfc, 



'-2k^M{k) 



2tt 
Qx)\{ky\Aj\qy)\'^Y 



(19) 



In the mean time, kpXc 3> 1 permits us to assume 
that the reflected/transmitted momenta q/q' at the inter- 
face are not much different from their values without the 
roughness. Then, only the q^ ~ kx survives in the brack- 
ets and can be pulled out of the summation to reproduce 
Eq.(3) in Rcf.[ll|. 



III. RESPONSE FUNCTION 

For a more systematic study of the IR in IXC, we de- 
note the Fourier component {ky\Aj\qy) by Hj[qy — ky) 
and all terms in Ea. ([T5)) and Eq. (fT6|) are proportional to 

Hj{qy ~ ky)Hl{ky - qy) . 



For simplicity, assume that the two side layers arc made 
of the same material and so the index j may be omitted 
in the scattering coefficients. Then, after changing the 
integration variable from qy into Py = qy — ky, we can 
rewrite Eg. ljTS)) as 
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and the momenta q^ and q[^ are functions of qy — Py + ky. 
The Py in G{Py) signifies the momentum transfer in- 
duced by a given Fourier component of the interface pro- 
file. The response function G{Py) is determined by the 
reflection coefficients, momenta and the exponential fac- 
tors from the round-trip propagations between interfaces. 
The study in the Casimir effect concluded [l^l that the 
perturbation method would become equivalent to PFA 
in the limit of long corrugation wavelengths. We shall 
prove in the following that this statement remains true 
for IXC. By use of the relation G{Py) = G{—Py) implied 
by Eq.dUl) and Eq.dlSl) and taking the limit Py -^ 0, the 
sum of Eq.dini) and Eq.(I21]) becomes 
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where the summation can be carried out to give 

^ iHLiPy) - HniPy) ' = ((Ai(y) - An{y)f). (26) 



Summarizing the above calculations, IR introduces a 
shift to the coupling energy AE'-^'> in Eq.([2]) for interfaces 
with mild corrugations: 
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under the hmit of A < 2Tr /kp <C Ac- Note that Eq. ip?)) 
is of the form of PFA which Taylor expands the varia- 
tion Aj^d/) — Aji{y) to the second order for the coupling 
energy. This demonstrates that Eq. ([20|) and Eq.([22|) 
are more general than PFA since they do not require 
kpXc » 1. In the special case of AL{y) = Aji{y) ^ 0, 
our method can still capture the effects of IR while PFA 
predicts none. 

In the next section, we apply our method to real 
systems where kpXc is not necessarily large, and con- 
centrate on the effect of correlation term Eg. ([20]) and 
Eg. ipT]) . Discrepancies between our results and those of 
PFA will be highlighted. The propagation term will be 
shown to display interesting features from the quantum 
interference for GMR. By increasing the potential barrier 
V{) of the spacer above Ep, our previous results can be 
applied to TMR. That turns the propagation term into 
a decaying function of the spacer width and the mathe- 
matical form of Eq. ([20l) through Eq. ([23| shall bear more 
resemblance to those of the Casimir effect llSl. Il9ll . 



IV. CORRELATED INTERFACES 

Let us start from simple sinusoidal functions for the 
corrugation on two correlated interfaces: 
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where p = 2tt/Xc. The energy correction thus depends 
on the lateral mismatch b. Plugging Eq. (p8)) into Eg. ([20)) 
gives 
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Similar to the procedures from Eg. ([25]) to Eg. ([27| . 
Eg. ip^ can be turned into the PFA form in the limit 
of p — > and 
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Equation (|29p implies the correlation energy correction 
can be modulated by a sinusoidal function of the phase 
difference between the two interfaces. Since the uncorre- 
lation term does not depend on the phase difference, the 
response function G'^{p) for the correlation term can be 
measured by substracting the coupling energy in Eq.([T3]) 
for in-phase case from that for out-of-phase. In the next 
two subsections, we shall use the sensitivity function 
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to quantify the discrepancy between our results and the 
PFA ones. 



A. TMR 

The sensitivity function p"^ for a typical TMR system 
is plotted in Fig[T] as a function of p/kp for different 
values of spacer width D. The p'^ in this figure exhibits 
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FIG. 1: (Color online) Sensitivity function p" is plotted as 
a function of p/kp for Vo/Ef = 2 and kpD — 3 (in dotted 
line), 7 (dashed line), and 10 (solid line). 

three traits: (1) it never exceeds unity which implies that 
the PFA always overestimates the correlation effect for 
TMR; (2) it approaches unity at small p/kp when our 
method reduces to the PFA; (3) it decays exponentially 
to zero when p/kp becomes large, which is corroborated 
by our analytic derivations for the asymptotic form of 
G'^{p) = ap e~P^ at p/kp ^^ 1 by Taylor expanding 
Eg. (pT|) . The parameter a depends on Ep and Vq. These 
features are shared by the Casimir effect [l^l because of 
their similar mathematical formalism. 



B. GMR 

As we reduce the potential barrier to Vq/Er — 0.5, the 
system enters the GMR regime. The results are plotted 
in Figl21 Comparing to FigH] for the TMR, the sensi- 
tivity function becomes oscillatory and can be negative 
in certain ranges of p. The period of oscillation shortens 
as the spacer gets thicker because the spacer width D is 
multiplied to the corrugation period p in the phase term. 
Furthermore, p'^ can now exceeds unity which is a neces- 
sary condition for the significant enhancement of IXC by 
IR[11|. Detail derivations for the enhancement and this 
statement are arranged in the following section. 

If the spacer potential is changed into a well, the 
emerging bound states is expected to have a limiting 
role at mediating IXC since their probability decay ex- 
ponentially into the metallic side layers. This is indeed 
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FIG. 2: (Color online) Variation of p'^ as a function of p/fc^ 
for VqJEf = 0.5 and kpD — 3 (in dotted line), 7 (dashed 
line), and 10 (solid line). 



true for smooth interfaces. However, in the presence of 
IR, the new eigenstates are a mixture of both scatter- 
ing and bound states. And, according to Eq. ([M|) . the 
bound states can contribute and render the denomina- 
tor, 1 — M{q), in Ea. (P5)) vanishing. So it is expected 
to generate new features in IXC. We divide the response 
function at Eq.([2T|) into two terms: 



'=• dE f dk 



GlJPy) = -21va -JL^k^q^ 



M{k)M{q) 



1 - M{k) 



(32) 



27r / 27r 



Gup.)^-.i..r'^i'^.kj''''''''^ 



'i2)K^y> --/ 27r7 27r 

1 M{q) + M{k) 

2 l-M{q) ' 



1 - M{k) 



(33) 



The first term is plotted in Fig|3{a), which represents the 
contribution from the paths that are scattered to momen- 
tum q in the beginning and then reflected to the original 
momentum k. The second term, as shown in Fig|3l^b) 
and is directly related to the additional resonance trans- 
mission, represents the contribution from the remaining 
paths that are transmitted in momentum g for more than 
one loop before being reflected to k. To clarify the dis- 
similarity of these two terms, we calculated a well system 
in FigU] and show the magnitudes of the first term, sec- 
ond term, and the total value of the response function. 
Again, the response function is divided by Gc{0) in the 
figure. Because the well system for a trilayer exhibits two 
bound states which generate two singularities in the in- 
tegral of Eq.(j331) for resonance transmissions, the second 
term of the response function displays two kinks as we 
modulate the corrugation period. 

Although we use a perturbative method for the calcu- 
lation of IXC energy between rough interfaces, the result 
of correction in GMR system still present unusual charac- 
teristics from the quantum interference and the resonance 
states. The PEA is not relevant to this regime. 



FIG. 3: (Color online) Diagrams for the scattering paths in 
(a) the first and (b) second terms of the response function in 
Ea. (|32[) and Ea.(|33[l. respectively. The thick/thin lines denote 
the paths with momenta q/k. (c) The dashed line denotes the 
path which is contributed by the thin line with any number 
of loops 
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FIG. 4: (Color online) Variation of p'^ versus b/kp with 
Vo/Ef = -0.2, kpD = 3 and for Eq.^ (dashed line), 
Ea. p3[) (thick line) and their sum (thin line). 



ENHANCEMENT OF IXC 



Equation (|13|) consists of two parts, correlation term 
in Eq.(15) and uncorrelated one in Eq.(16). By use of 
Eq. (f27|) from PEA, IXC for mildly corrugated interfaces 
can be expressed as: 
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where the Aj^/ji and D arc the same definitions as in 
previous sections. Since AE'-'^'f for GMR is an oscillatory 
function of the spacer width D, it shares the same sign 
as the negative of its second derivative for kpD > 1. As 
a result, the second term always strengthens IXC, while 
the third term diminishes it. An overall enhancement 
of the coupling strength is realized when the correlation 
term dominates. However, this is not possible in the 



above PFA expression for mild corrugations. The sum of 
these two terms can never be positive and, at most, cancel 
each other to give null contribution when the tomography 
on both interfaces happen to be unrealistically identical. 
Therefore, it is safe to say that the IR also suppresses 
IXC within the second-order perturbation of PFA. 

We shall now demonstrate that more severe corruga- 
tions, A ^ 2Tr/kF ~ Ac, and correlated tomographies 
are two essential ingredients to enhance IXC. The former 
requires us to improve upon the PFA within the second- 
order perturbation, while the latter brings in strong 
quantum interference. To clarify this statement, we use 
our method to estimate the coupling strength for the 2D 
trilayer system with Vo/Ep = —0.2 and identical topog- 
raphy A{y) = l/(2fcF) sin(fcFy/2) on both interfaces. As 
shown in FiglSl the effect of correlated and severe cor- 
rugations can improve IXC for smooth interfaces by as 
much as one and a half times. 
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FIG. 5: (Color online) Coupling strength is plotted as a func- 
tion of kpD for two interfaces that share the same topography 
A{y) = l/(2fcF) sin(fcFy/2) with Vq/Ef = -0.2. The result 
for smooth interfaces is shown in the dashed line for compar- 
ison. 

It is heuristic to approximate (PAE^^^/dD^ by 
-4,kpE^°'> when kpD > 1 since IXC for two identical 
interfaces with topography A(y) can be neatly reduced 
to 
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It is then clear that an enhancement in IXC is being 
caused by the dominance of the correlation sensitivity 
function over the uncorrelation one, which can be real- 
ized for a wide range of p/kp in FiglHl The largest en- 
hancement appears around p = 0.22fci?, which is about 
two times that for smooth interfaces. We checked that 
Eq. ([55)1 gave roughly the same value as that without the 
approximation, 1.4 and 1.5 respectively for the topogra- 
phy and parameter in Figl5l 

We have demonstrated that our approach is more gen- 
eral than the PFA at capturing the effect of quantum 
interference among different reflected paths of carriers 
within the spacer. This allows us to obtain the enhance- 
ment of IXC when p is comparable to kp. It is then 
reasonable to ask what happens when the characteris- 
tic corrugation p is much greater than kp. It turns out 




FIG. 6: (Color online) Setting the corrugation function A(y) 
to observe the sinusoidal form, sinpy, the correlation p'^ (solid 
line) and uncorrelation sensitivity functions p"'^ (dashed line) 
are plotted as a function of p/kp for Vo/ Ep = —0.2 and 
kpD = 2.2. 



that the momentum change along the surface of the in- 
terfaces, being of order p, is so large that the longitu- 
dinal momentum becomes pure imaginary during elastic 
collisions. This means that these paths only survive a 
short distance in the a;-direction and, therefore, are not 
expected to lead to major interference. To be more rigor- 
ous, if we approximate the additional scattering momenta 
Qx, q'x due to IR by ip and insert into Eq.(21) and Eq.(23), 
the sensitivity function for correlation term will decay as 
p e~P^ while the uncorrelation one remains roughly a 
constant. In retrospect, the process of allowing different 
reflected paths of carriers to interfere is similar to that 
of localization, but the eventual effect is different. In 
this case of IXC, although the second term in Eq. p5|) 
contains the second power of the small perturbation pa- 
rameter kpA, the quantum interference is still capable of 
rendering this term large via the other factor, p'^ — p^'^. 



VI. DISCUSSIONS AND CONCLUSIONS 

Our perturbative approach to evaluate the effects of 
interface roughness in the trilayers was motivated by a 
similar effort in the Casimir problem. To be precise, the 
role of interface roughness at causing interference among 
reflected electromagnetic waves within the cavity finds a 
better analogy in TMR, rather than GMR. Reason be- 
ing that the decaying nature of carrier wavefunctions in 
TMR limits their quantum interference within the spacer, 
while the fact that virtual photons lack a characteristic 
length scale like the inverse of Fermi momentum l/kp 
dilutes any possible constructive interference after all 
wavelengths are summed over. In contrast, the quan- 
tum interference survives and has a nonnegligible effect 
on GMR when the characteristic wavelength of corru- 
gations is shared by both interfaces and comparable to 
l/kp. 

Compared to our previous study [ll|, a couple of im- 
provements have been made in this work. First, although 
both calculations retained up to the second order in the 



corrugation amplitude A(y), we included more loops of 
multiple scattering from the smooth part of the inter- 
faces. Furthermore, A{y) was no longer confined to be of 
the sinusoidal form. Second, an analogy to the Casimir 
effect was made, which allowed us to borrow the con- 
cept of response and sensitivity functions as indicators of 
the extent of influence by A{y) without having to know 
its detailed form. We followed up by more discussions 
on similar and different effects of A(y) in trilaycrs and 
the Casimir mirrors. Third, more thorough derivations 
were done to compare our approach with the prevail- 
ing proximity-force approximation for different periods 
of corrugations. This enabled a better quantitative esti- 
mate of the enhancement from quantum interference. 

In conclusion, we find that the perturbative approach 
reaches the same conclusions as the proximity-force ap- 
proximation in the limit oi p ^ kp. Namely, mild cor- 
rugations lead to a suppression of the interference and 
thus the interlayer exchange coupling. Correlated rough- 
ness with short wavelengths gives rise to several inter- 
esting features: (1) The energy correction oscillates as 
we vary the corrugation wavelength. (2) The magnitude 
of correction can be larger than the prediction made by 
the proximity-force approximation. (3) While they are 
expected by the proximity-force approximation to be ir- 
relevant to the transmission coefficient, the bound states 
within the spacers are found to affect the resonance trans- 
mission through several kinks in the energy correction. 
One last important feature concerns the enhancement 
of interlayer exchange coupling by the interface rough- 
ness. Its occurrence relies on further requirement that 
the Fourier conjugates alluded to above be close to the 
Fermi momentum. 

Support by the National Science Council in Taiwan un- 
der Grant No. 98-2112-M007-005-MY3 is acknowledged. 



Appendix 

Appendix A: Connection between The Casimir 
effect and IXC 

In this appendix, we would like to extend the concept 
of radiation pressure in the Casimir effect to the IXC 
problem. We start from the Casimir energy which is the 
summation of zero-point energies for quantum states in 
the presence of boundaries: 



^= oll^" 



(Al) 



where N is total number of carriers and Er is the refer- 
ence energy without the boundary. 



1. Radiation Force in One Dimension 

To clarity the connection between these two energies, 
we now start from the radiation force of fields in a one- 
dimensional potential well with width D. The cigenstatcs 
of carriers in a quantum well can be written as a combi- 
nation of two travelling waves in opposite directions: 



— sm fc„a: = —i\ 
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(A3) 



where fc„ = rnr/D is the quantized momentum. 

The radiation force measures the impetus per unit time 
contributed by the carriers in the cavity: 
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= E^(^c-^L + ^c-^L) (A4) 



where w„ — kn/m and m denotes the carrier mass. The 
IXC energy can be evaluated by 



AE^- dD'AF{D') 

J oo 

N 

= ^En{^,G-etil; + Kc-Kc)~Er- (A5) 



It is not surprising that the above expression reduces to 
Eq.(|A2| upon assigning (etic ■ etc) ^ (^c ■ ^c) = 1/2. 
Similar derivations for the vaccuum states give: 
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and 



where n and w„ denote the n-th bound state with fre- 
quency LOn- In contrast, the IXC energy measures the 
increase in the total carrier energy when considering the 
boundaries: 
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E = ^i^n{et.c -e^c + ^c ' ^c)vac = 9 E'^"- (A'^) 



where the extra coefficient 1/2 comes from the fact that 
the quantum amplitude for the vaccuum state is just 
one half of the corresponding commutator from photon 
operators J25J. 
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Based on the above discussions and the relations T< = 
T>, |-R>P + |2^>P = 1, the quantum radiation force in 



Eq. ()A8|) can be rearranged as 
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D^\^)-2 (A12) 



where jC^p and \D^\'^ are functions of {Tl,Tb.,Rl,Rr} 
which denote the transmission and reflection coefficients 
at each of the barriers in FiglT] Equation (jA12[) will 
become 



FIG. 7: Schematic plot of a trilayer system to clarify our 
notations: L/R denote the fields in the left/right sides of the 
spacer, whereas C represents those in the spacer. 



Radiation Force in Three Dimensions 



We generalize the result in Eq. (jA4[) to estimate the 
radiation force in a 3D trilayer system in FiglT] 
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where L/R denote the fields in the left/right sides of 
the trilayer and L^ is the system length in x direction. 
The extra factor cos^ comes from projections of the 
momenta and velocities on the normal direction of inter- 
faces. Take the continuum limit and the summation can 
be changed to an integral 
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where A = LyL^ and the notation IBZ signalizes the 
range of integral to be bounded by the interfacial Bril- 
louin zone. The amplitude of quantum fields can be ob- 
tained from the scattering states in the trilayer system: 



i^k,>(r) 



where > denotes the moving direction of the scattering 
state. The left-moving one ipk,< can be defined similarly. 
Since these two states are orthogonal, their contributions 
to the inner product of field amplitudes can be separated: 
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We can also calculate the IXC energy: 
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Equation (jA14p . which has been proved in Ref.[a] by 
using the concept of quantum interference and Green's 
function, is commonly used in the study of IXC. 

In the Casimir problem, we follow the same procedures 
for the vaccuum state in Eq.(|A7| to obtain 
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where p denotes the transverse electric and magnetic 
modes. By use of the Cauchy theorem, we can shift the 
integration to the imaginary frequency axis and rewrite 
the Casimir force in Eq. (JAISP and energy in Eq. (JA16I) as 
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where k = ./fcjf +S.'^/c^. Equation (|A17p and Eq. (|A18|) 
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(All) 



1 

are the familiar formulae for lossy optical cavities. A 
previous article |25| has used the vaccuum radiation pres- 
sure to derive Eq. (JA17I) . Although similar in concepts to 



theirs, our derivations to relate the IXC and the Casimir 
energy are more straightforward. One quahtative differ- 
ence is that the upper bound of energy integration in 



IXC is bound by the Fermi energy which energy scale 
eventually renders the IXC force oscillatory in 2kpD. 
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